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Abstract 

The Menter Shear Stress Transport (SST) k — uj turbulence model is one of the most 
widely used two-equation Reynolds-averaged Navier-Stokes turbulence models for aerody- 
namic analyses. The model extends Menter’s baseline (BSL) model to include a limiter 
that prevents the calculated turbulent shear stress from exceeding a prescribed fraction 
of the turbulent kinetic energy via a proportionality constant, ai, set to 0.31. Compared 
to other turbulence models, the SST model yields superior predictions of mild adverse 
pressure gradient flows including those with small separations. In shock - boundary layer 
interaction regions, the SST model produces separations that are too large while the BSL 
model is on the other extreme, predicting separations that are too small. In this paper, 
changing a± to a value near 0.355 is shown to significantly improve predictions of shock 
separated flows. Several cases are examined computationally and experimental data is also 
considered to justify raising the value of a\ used for shock separated flows. 


Nomenclature 


Cf 

Clim 

d 

H 

k 

Vk 

Re x 

Sij 

t 

Uj 

Coo 

-?/+ 


Wij 

y + 

e 

K 


P 

Pt 

UJ 

n 

p 


skin friction coefficient 

shear stress limiting coefficient 

distance from wall 

backward facing step height 

turbulent kinetic energy 

production of turbulent kinetic energy 

plate Reynolds number 

mean rate of strain tensor 

time 

mean velocity vector 

mean freestream velocity 

velocity normalized with wall shear stress 

fluctuating velocity vector 

fluctuating velocity components 

mean vorticity tensor 

Cartesian coordinates 

wall normal coordinate 

boundary layer thickness for UFAST test case just upstream of shock generator 

turbulent dissipation rate 

von Karman constant 

dynamic viscosity 

dynamic eddy viscosity 

specific turbulent dissipation rate = e/k 

vorticity magnitude 

density 

turbulent stress tensor 


NAS A/TM— 20 13-217851 


1 



Introduction 


Despite extensive efforts in the aerodynamics community over more than two decades, the prediction of 
flow separation remains elusive. In supersonic flows, the shock wave / turbulent boundary layer interaction 
(SWTBLI) produces especially complex flow separations that have significant effects on the operability 
of specific components such as aircraft inlets and also on overall performance of aerospace vehicles. A 
recent workshop considering CFD calculations for a set of SWTBLI cases was held in conjunction with 
the 48th American Institute of Aeronautics and Astronautics (AIAA) Aerospace Sciences Meeting. The 
focus of this workshop 1,2 was SWTBLI as occurs in supersonic inlets. Several investigators contributed 
solutions to this workshop including both Reynolds-averaged Navier- Stokes (RANS) computations and Large- 
Eddy Simulations (LES) computations. DeBonis et al 3 provides a comprehensive assessment of the CFD 
calculations, including uncertainty analysis of the submitted CFD results and experimental data obtained 
for the same configurations. The results of the workshop indicated that none of the RANS methods were 
able to capture the details of the SWTBLIs and that LES, while offering promise for the future, still cannot 
be applied to configurations at realistic Reynolds numbers. 

Of the RANS turbulence models applied to the SWTBLI cases in the workshop, the most commonly 
used model by the participants was the two-equation Shear Stress Transport (SST) model of Menter. 4 The 
SST model is a k-u; formulation that has become one of the most widely used turbulence models for RANS 
analyses. It has been shown to be robust and relatively accurate for a broad range of flows, including 
wall boundary layers and free shear layers. For small separations due to mild adverse pressure gradients 
in subsonic flows, the model has demonstrated success. As observed in the results of the workshop and 
elsewhere, however, the SST model tends to over-predict the size of flow separations in SWTBLI problems. 
In Ref. 4, Menter also introduced a similar k-cj model without a stress limiter, called the baseline (BSL) 
model, which conversely tends to underpredict the magnitude of shock-induced separations. 

Because the BSL and SST formulations predict significant differences in shock induced flow separations, 
on opposite extremes of the experimental data, we investigated the distinguishing aspect of the SST model, 
the shear stress limiter, in this paper. Details of the Menter SST model and background of the shear stress 
limiter that has made the model so popular are described. It is shown that for flows with shock-induced 
separations, a minor adjustment to the key coefficient in the shear stress limiter enables improved predictions 
of flows containing SWTBLIs. 

Menter SST Model 

For nearly two decades, the practical state-of-the-art in RANS turbulence modeling has been two-equation 
eddy viscosity turbulence models. In these linear eddy viscosity models, the Boussinesq approximation is 
used to relate the turbulent Reynolds stress to the mean rate of strain tensor through a turbulent (or eddy) 
viscosity. The SST model is a two-layer eddy viscosity model which employs the original (commonly referred 
to as the 1988 version) of the k — uj model of Wilcox 5 in the inner region of boundary layers and switches 
to a k — e model in the outer region of boundary layers and in mixing regions. The outer k — e model is 
transformed to the k — uj form with a blending function used to transition between the two layers. 
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The constants for the inner mode (set 1) are as follows: 


(10) 

( 11 ) 


71 = ^ - ““TpT’ £7*1=0.85, o-a,! =0.5, ft = 0.075, 

and the constants ( 02 ) for the transformed k — e model (set 2) are: 

72 = g-^, <7*2 = 1.00, a w2 = 0.856, ft = 0.0828, 

ft Vft 

These constants are blended using the same switching function, F\ that is found in the model equations 
such that 0 = Fi0i + (1 — T\)02 for any of the given parameters. The remaining constants have the following 
values: 


0* = 0.09, k = 0.41, ai = 0.31, 

The SST turbulence model has become popular for its accuracy and applicability to complex geometry 
due to the strengths of the k — uj model in the near wall region, and k — e model in free shear layers. The 
BSL model also introduced by Menter, 4 is identical to the SST model with two exceptions. First, the value 
for Oki in set 1 is set to 0.5, which is the same as in the original Wilcox model. 5 Secondly, the turbulent 
viscosity is simply /q = pk/co. Referring back to Eqn. 9, the SST model employs a “shear stress limiter” to 
prevent the calculated turbulent shear stress from exceeding a prescribed fraction of the turbulent kinetic 
energy. It is this term that we focus on in this paper, for while the SST model has demonstrated success 
for mild adverse pressure gradient flows, the prediction of SWTBLIs has generally resulted in shock-induced 
separations that are much larger than those measured in experiments. 

Shear Stress Limiter 

Menter 4 remarked that one of the primary deficiencies of eddy- viscosity models, in comparison to full 
Reynolds- stress closures, is that they do not account for the transport of the turbulent shear stress. Menter 
refers to the “Bradshaw Assumption” where the turbulent shear stress in a boundary layer is proportional 
to the turbulent kinetic energy through a constant as shown in Eqn. 12: 
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r = — pu'v ' = paik (12) 

The proportionality constant, ai, used in the limiter for the SST model, is also commonly referred to as 
the structure parameter, or the “Townsend 6 structure parameter.” In the literature, the former is usually 
defined as — u'v' = a\q where q = 2 &, and a\ will be half of the value corresponding to that used in Eqn. 12. 
Here, we will refer to the form denoted by Eqn. 12 when discussing the limiter or the experimental findings 
in the following sections. Because two-equation turbulence models calculate the turbulent shear stress from 
Eqn. 4 or the nearly approximate form using the vorticity magnitude, U, Eqn. 9 limits the eddy viscosity 
calculated by the model to the smaller of pk/uj and pa\k/VtF 2 . As shown in Ref. 7 the F 2 function is unity 
in the inner three-fourths of a zero-pressure gradient boundary layer and drops to zero outside boundary 
layers so that no limiting is performed outside boundary layers. 

It is obvious then that the selection for the value of the proportionality constant, ai, is crucial to the 
performance of the SST model. Menter selected a\ = 0.31 as the standard value for the SST model. Coakley 8 
developed a q — uo model which also used a shear stress limiter with a value for this constant of 0.30. Both 
the Coakley and Menter models followed the work of Bradshaw et al. 9 where a one-equation model was 
formulated based on the same relation of the turbulent shear stress to the turbulent kinetic energy via 
a structure parameter set to approximately 0.30. Bradshaw et al. 9 and Bradshaw 10 utilize experimental 
data of Klebanoff 11 to illustrate that this value of the structure parameter is characteristic of zero pressure 
gradient flat plate boundary layers. For mild adverse pressure gradients, Aubertine and Eaton 12 as well as 
Alving and Fernholz 13 show that as the flow passes through an adverse pressure gradient and even small 
separation, the structure parameter actually drops from the zero-pressure gradient value, before recovering 
further downstream of reattachment. 

The success of the SST model using a\ = 0.31 relative to the otherwise similar BSL model, has been 
well demonstrated by Menter 4 and many other investigators for attached boundary layers or mild adverse 
pressure gradient flows leading to relatively small separations. The motivation here is to examine the effect 
of the shear stress limiter, and the original value of a\ that was set to 0.31, for flows outside these regimes 
and in particular to focus on the SWTBLI problem. 

The value of a\ for a homogenous free shear layer flow investigated by Rose 14 was found to be approx- 
imately 0.35. Bradshaw 10 concluded that the empirical relationship between turbulent shear stress and 
turbulent kinetic energy for boundary layers that are attached, or nearly so, is not likely valid for other 
turbulent flows. The fact that the SST model employs the F \ function enables the shear stress limiter to be 
turned off outside boundary layers, which is an advantage of the hybrid model formulation. 

Other reported efforts have investigated a different value than a\ = 0.31 for shock-boundary layer dom- 
inated flows. In the most recent version of Wilcox’s k — uo model , 15,16 commonly referred to as the Wilcox 
2006 k — uo model, a shear stress limiter is introduced. This limiter differs from Menter’s formulation in that 
it uses (1) the rate-of-strain tensor instead of vorticity, and (2) a coefficient Cum rather than a\. For shear 
dominated flows, the effective structure parameter is given by Eqn. 13. 



Using Wilcox’s model coefficients yields a value of a\ = 0.342. Wilcox also mentions that he examined 
using the vorticity magnitude with effective a\ set to about 0.316. 

Another significant difference between the Wilcox and Menter models is that this limiter is applied 
everywhere, while in the Menter SST model the F 2 function only applies the shear stress limiter in the 
boundary layer. Tan and Jin 17 investigated the SST model with variations in a\ also expressed through the 
context of the Cu m parameter and concluded that a constant value for a\ was not likely. Edwards et al. 18 
investigated variations in the effective a\ used by SST as the inner model portion of a hybrid RANS-LES 
formulation and concluded that a\ = 0.356 provided the best agreement with experimental results for a 
Mach 5 compression corner problem. 

It should also be noted that algebraic Reynolds stress models inherently act to reduce the Reynolds 
shear stress in similar circumstances. This stress limiter is not a hard limit like those discussed above, but 
is instead a function of local flow invariants and the coefficients used to model higher-order terms such as 
the pressure-strain correlation. Wilcox 16 provides a brief analysis of results from this class of model. Such 
modeling efforts are beyond the scope of the present work, but may hold promise for the future. 
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While the work reported in this paper is not the first study to investigate the effect of the shear stress 
limiter on calculation of shock separated flows, it is the first study we are aware of that examines the details 
of turbulent statistics measured in the experiments alongside the computational results to determine a more 
appropriate setting for the stress limiter constant in flows where shock-boundary layer interactions are a key 
flow feature. 

Structure Parameter in SWTBLI Regions 

For SWTBLI flows, there is experimental evidence that the structure parameter can significantly exceed the 
values that are typical of the aforementioned boundary layer problems. Smits and Muck 19 present results 
for a series of compression corner flows and note that the peak value for the structure parameter a " given 
by —u'v'/u' 2 may increase by over a factor of two with the peak location occurring in the outer part of the 
boundary layer. Smits and Muck, 19 Jayaram et ah, 20 and Kuntz et al. 21 show that the SWTBLI creates 
a large amplification in the turbulent stresses, which increases with compression corner turning angle for a 
fixed incoming freestream Mach number. 

In our experience applying the Menter BSL and SST formulations to SWTBLI problems of the 2010 
AIAA workshop 22 as well as that of other participating groups summarized in Ref. 3, it was generally found 
that the BSL and SST models provided predictions on either side of the experimental measurements. In 
particular, the BSL model underpredicted the magnitude of the separation while the SST model, using the 
standard a\ value of 0.31 significantly overpredicted the size of the separation. The objective of the present 
work is to determine if there is a more optimal choice for the value of a\ used within the context of the 
Menter SST model. Clearly the success of this model is largely due to use of the shear stress limiter for many 
flows. However, as illustrated by the experimental observations discussed in this section, there is evidence 
that limiting a\ to 0.31 for SWTBLI dominated flows may prohibit correct values for the turbulent shear 
stress to be modeled. 


Results 

In this section, we first examine variations in the a\ parameter for an incompressible flat plate flow, low 
speed axisymmetric diffuser, and incompressible backward facing step. Focus is then moved to SWTBLI, 
which is the target flow problem of interest. In the process of examining these flows, it was found that 
values of a\ greater than 0.40 were very similar to the unbounded BSL model. Therefore, the range of a\ 
investigated here was restricted to 0.31 < a\ < 0.40. 

Flat Plate Boundary Layer 

A nearly incompressible flat plate flow was investigated with freestream Mach number set to 0.2 in order to 
provide reasonable convergence with the compressible flow solver. For this case and the next two benchmark 
low subsonic test cases, grids and boundary conditions that are available on the turbulence modeling website 
described in Ref. 23 were employed. All of the grids were previously verified to yield grid independent 
solutions. 

In addition to variations in ai, it was noted in an earlier section that cr^i for the inner model also differs 
between the BSL and SST models. For completeness then, a BSL case was run to investigate the effect of 
this coefficient, which controls the diffusion of turbulent kinetic energy. The skin friction coefficient along 
the flat plate is shown in Fig. 1 and velocity profile is shown in Fig. 2. Both are compared with data of 
Wieghardt and Tillman. 24 The line denoted by only “SST” in the legend of this and subsequent figures has 
a\ = 0.31 as in the default model. The skin friction plot shows that there is minimal change in the BSL 
results due to changing (Jk i from 0.5 to 0.85. 

More details of the actual turbulence model behavior are shown in the plots of eddy viscosity in Fig. 
3 and turbulent shear stress normalized by the turbulent kinetic energy (or effective structure parameter, 
ai) in Fig. 4. The experimental data in Fig. 4 is taken from Klebanoff 11 but is for nearly an identical 
plate Reynolds number, 4.2e6, as for the mean velocity profile in Fig. 2. The SST model produces the 
smallest values of eddy viscosity, which is due to the limiting behavior when using a\ = 0.31. The inflection 
in the eddy viscosity and discontinuity in the slope of the turbulent shear stress for the SST solution near 
y+ = 2000 is due to the switch from inner model to outer model via the F 2 function. The “SST” solutions 
where a\ is 0.34 or higher are indistinguishable because the shear stress limiter is never reached. These 
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solutions also match the BSL solution with <Jk 1 changed to 0.85, as expected. The BSL solution with <Jki 
set to its standard value of 0.5 allows for less diffusion of turbulent kinetic energy, also as expected. 

Axisymmetric Separated Diffuser 

Perhaps the case which most clearly made a distinction between the SST formulation and other two-equation 
models in Menter 4 was the Driver axisymmetric diffuser. 25,26 In particular, the advantage of using a shear 
stress limiter in conjunction with the Menter SST model was demonstrated. This case considered a Mach 
0.088 flow over a cylinder aligned with the flow, with an adverse pressure gradient established by diverging 
the tunnel walls. The case is typically run in axisymmetric mode by establishing an outer streamline believed 
to accurately replicate the otherwise three-dimensional flow case. 

Utilizing the verified grid available from Rumsey et ah, 23 calculations were performed using the original 
SST and BSL formulations as well as SST with variations in a\. The pressure distribution plot in Fig. 5 
shows nearly identical results to those originally shown by Menter 4 and clearly shows the superior result 
enabled by the SST solution with the standard value for a\ compared to all of the other solutions which 
nearly collapse to a single line. The skin friction comparison in Fig. 6 also shows that the Menter SST model 
case differs from the others, which are all collapsing to the same solution, and indicates the largest extent of 
flow separation. 

With extensive turbulence data available from the Driver experiments, the details of the turbulent flow 
development through the adverse pressure gradient region were investigated as shown in Figs. 7-9. The 
turbulent shear stress contours in Fig. 7 show that the SST solution delays the growth of — u'v' compared 
to the BSL solution and the solutions with varying a\, particularly in the region 0 mm < x < 100 mm. This 
is the same region that, in Fig. 5, all of the pressure distributions from solutions other than the SST model 
begin to depart from the experimental data. One might conclude that a\ = 0.31 utilized by the original 
SST model is capturing the flow physics most accurately. However, by examining contours of —u'v'/k 
(Fig. 8), which is equal to one can clearly see that none of the solutions, including the original SST 
model, accurately reproduce the experimental measurements. The experimental results show that —u'v'/k 
drops slightly near the beginning of the adverse pressure gradient region ( x = —300 mm) and begins to 
increase near the end of the pressure gradient region ( x = 200 mm) . This trend matches that found in other 
experiments. 12, 13 The Driver data, however, indicate that —u'v'/k remains substantially below 0.31 until 
well beyond the end of the domain. 

The SST solution indeed limits —u'v'/k to 0.31 through most of the boundary layer, with a streak 
of —u'v'/k exceeding 0.31 where the model switches to the outer formulation and does not have a limit 
on —u'v'/k. The SST solutions with a\ > 0.31 and BSL solution allow for higher —u'v'/k lower in the 
boundary layer, as expected. The experimental measurements of turbulent kinetic energy in Fig. 9 show a 
much earlier increase over upstream levels (near x = 100 mm) than any of the computational results. The 
turbulent shear stresses utilized in the original SST model that enable the best pressure distribution, result 
from applying a fortuitous limit (0.31) on a turbulent kinetic energy field that has responded too slowly to 
the adverse pressure gradient in comparison to the experimental data. 

Backward Facing Step 

We next examine the backward facing step case of Driver and Seegmiller. 27 This case had flow at Mach 
0.128 approaching the backward facing step region. Once again a grid taken from the turbulence modeling 
resource website was utilized here. A comparison of experimental measurements and computational results 
for skin friction in the vicinity of the step is shown in Fig. 10. There is little difference between the solutions 
when ai is increased to 0.34 and larger values. The slight difference in the BSL solution from the SST 
solutions with increased values for a\ is due to the different value used for than any of the SST-based 
solutions. Eddy viscosity contours for the two extreme cases, SST and BSL, are provided in Fig. 11. 

In the experiments, laser Doppler velocimetry was used to obtain extensive data including turbulent 
stresses. Using their data which included u',v ', and —u'v', an estimate for the turbulent kinetic energy can 
be made by assuming that V w' 2 = \ ^v^F 2 + Vu 2 ^. Figure 12 shows the calculated effective structure 

parameter (—u'v'/k) from the calculations and the data of Driver and Seegmiller. The experimental data 
shows a streak of —u'v'/k of approximately 0.35 extends across the shear line originating from the step to 
x/H = 2.5, and while the peak value drops going further downstream, a streak where —u'v'/k exceeds 0.31 
extends to x/H = 8. 
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Focusing on the region in the separation beyond the step, one may observe that the SST solution with 
a\ = 0.355 and BSL solution predict slightly higher levels of —u'v'/k near 2 < x/H < 4 than the standard 
SST solution. However, none agree well with the experimental data which indicate higher levels of —u'v'/k 
near the wall. This likely leads to the differences between experiment and computation for skin friction 
coefficient in the separated region. 

Mach 2.25 SWTBLI 

As stated in the introduction, the primary problem for which we suspect that limiting cq has had a detri- 
mental effect is flows involving SWTBLIs. One of two SWTBLI configurations considered in the previously 
mentioned AIAA workshop was a test case obtained at the Institut Universitaire des Systemes Thermiques 
Industriels (IUSTI) in Marseille, France. 28 The experimental data was part of the European Union SWT- 
BLI research project referred to as UFAST. 29 The UFAST experiments utilized an 8 degree shock generator 
which spanned the entire width of the tunnel with an approaching Mach 2.25 flow. The supply stagnation 
pressure was 50.5 kPa and the stagnation temperature was 293 K. Particle Image Velocimetry (PIV) was 
utilized to measure mean velocities and also obtain turbulent statistics. A three-dimensional computational 
grid was constructed having approximately 7.4 million points, and the average y+ of the first point off of 
the wall set to 1. A grid sensitivity study indicated minimal differences between solutions obtained with this 
grid and solutions using a sequenced grid, that is using every other point in each of the three computational 
directions. 

In addition to simulations with the Menter BSL and SST models (using a\ = 0.31), several other values 
for a\ were investigated. A comparison of axial velocity contours from these simulations and experimental 
measurements is shown in Fig. 13. It may be observed that in all cases, the SWTBLI is centered around 
x — 320mm. Clearly the unmodified SST model provides an interaction region that is much larger than the 
experimental data and the BSL model indicates an interaction region that is too small. Table 1 provides a 
measure of the separation extent for the experimental data and CFD solutions that is obtained by examining 
where the axial velocity is negative at two vertical positions from the wall: y\ = 0.187mm = 0.017£ w and 
i )2 = 0.813mm = 0.074£ u . Examining these results and the contours in Fig. 13, it appears the solution 
obtained for a\ = 0.355 provides the closest agreement with experimental data. Further comparison of axial 
velocity profiles in Fig. 14 supports these findings but also shows that none of the solutions agree well with 
the experimental data toward the end of the SWTBLI, where the flow is recovering to an attached state. 

Table 1. Length of negative axial velocity streak for UFAST SWTBLI case. The vertical positions where the 
separation lengths were measured correspond to y\ = 0.187 mm = 0.017(5 n and 1/2 = 0.813 mm = 0.074 8 U * 


Case 

Separation Length- 1 (mm) 

Separation Length-2 (mm) 

Experimental Data 

22.3 

15.1 

Menter SST, a x = 0.31 

52.9 

43.1 

Menter SST, a x = 0.34 

30.6 

19.6 

Menter SST, a\ = 0.355 

24.6 

13.5 

Menter SST, a\ = 0.37 

21.2 

9.7 

Menter BSL 

16.0 

0 


The experimental data suggests that raising a\ closer to 0.355 is warranted. Figure 15 shows the behavior 
of —u'v'/ k in the SWTBLI using the UFAST experimental data. As was the case for the backward facing step, 

V w' 2 was approximated as the mean of \fu and a/U 2 . Upstream of the interaction region, the attached 
boundary layer reaches a peak value for —u'v'/k that is slightly less than 0.30, which is in agreement with 
the boundary layer findings mentioned earlier. In both the initial portion of the SWTBLI and then towards 
the latter half of the separation region x > 320mm, it may be observed that —u'v'/k increases well beyond 
the values characteristic of equilibrium boundary layers. It is interesting to note that for the profile nearest 
the center of the interaction, the peak value of —u'v'/k temporarily drops very low. It is not surprising then 
that the unmodified SST model with a\ = 0.31 limits the turbulent shear stresses too much and results in 
a separated flow region that is too large. Humble et al. 30 took PIV measurements for a similar SWTBLI 
flow problem and also observed similar increases in the turbulent shear stresses in the latter half of the 
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SWTBLI and further downstream. These finding also agree with those of Refs. 19-21 which concluded that 
the turbulent stresses are amplified in the outer portion of the recovering boundary layer in SWTBLIs. 

Mach 5 SWTBLI 

We next examine the Mach 5 SWTBLI investigated experimentally by Schulein. 31 Due to the significantly 
higher Mach number of this case, a detailed mapping of turbulent statistics was not available. Comparisons 
are limited to the skin friction through the SWTBLI. 

The effect of variations in a\ on the computed skin friction is provided in Fig. 16. None of the solutions 
match the experimental data very closely. All of the solutions predict a more rapid drop in skin friction 
coefficient than the experimental measurements in the region where the oblique shock intersects the boundary 
layer. Similar to the previously considered Mach 2.25 SWTBLI case, it is obvious again that the SST solution 
indicates the largest separation region and the BSL provides the smallest, with the same trends in separation 
extent resulting from variation in a\. 

Scramjet Internal Flowpath 

The last case considered in this paper is the internal flowpath of a ground test scramjet configuration known 
as the Hypersonic International Flight Research Experimentation (HIFiRE) Direct-Connect Rig (HDCR). 
The configuration is a dual-mode scramjet test article with a long isolator, cavity flame holder, and multiple 
sets of fuel injectors. The case presented here simulates a flight Mach number of 5.8 and has an overall fuel 
equivalence ratio of 0.65. Due to the harsh environment of this reacting flow configuration, only centerline 
pressure distributions are available from Ref. 32 for comparison. 

A set of reacting flow test cases are compared. There are many factors in this case beyond SWTBLI that 
affect the resulting pressure distributions. These include thermal boundary layer effects, turbulent-chemistry 
interactions, and compressibility. However, as may be observed in Fig. 17(a), there is a significant SWTBLI 
region at the end of the isolator 8m < x < 15m which is where attention is focused here. The core flow 
Mach number approaching this region is 2.5. The comparison of pressure distributions in Fig. 17(b) shows 
a large discrepancy between the BSL and SST limits. For clarity, not all of the solutions are shown in the 
plot, but the trend is increasing pressure through the SWTBLI at the end of the isolator and beginning 
of the combustor as ai is increased, and the turbulent shear stress is permitted to be higher. The best 
agreement was provided by the case with a\ = 0.355. While not shown here, a second case was examined for 
a freestream Mach number of 8.0, equivalence ratio = 1.0, and core flow Mach number = 3.5 approaching 
the end of the isolator. The behavior of the turbulence models was very similar with a\ — 0.355 providing 
the best agreement. 


Conclusions 

The Menter Shear Stress Transport (SST) k — uj turbulence model has been one of the most successful and 
popular Reynolds- averaged Navier-Stokes turbulence models during the past two decades. A large reason for 
its success is in predicting wall bounded flows for a variety of flow conditions. Central to the model is the use 
of a limiter, which in the standard published version of the model prevents the turbulent shear stress from 
exceeding 0.31 times the turbulent kinetic energy in the inner part of a boundary layer. Menter’s baseline 
(BSL) model does not limit the shear stress. Experience with the two otherwise nearly identical models has 
shown that for flows dominated by shock wave / turbulent boundary layer interactions (SWTBLIs), the SST 
model typically overestimates the size of flow separations while the BSL model underestimates separation 
magnitude. In this work, we examined the effect of the proportionality constant, a\ which is used by the 
SST model to limit the turbulent shear stress. 

A set of several nearly incompressible cases were investigated first, including a flat plate, axisymmetric 
diffuser, and backward facing step. In addition to examining the computational solutions, the details of the 
experimental results were also examined with focus placed on the turbulent shear stresses, turbulent kinetic 
energy, and their ratio, which is referred to as the structure parameter in the literature, and equivalent to a\. 
For the low speed separating axisymmetric diffuser flow, the original SST model with a\ = 0.31 yielded the 
best agreement in terms of pressure coefficient through the separation region. However, the experimental 
turbulence measurements for this case indicated that a\ was substantially below this limiting value through 
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the separation, and that the best agreement enabled by the SST model with a\ = 0.31 may be considered 
fortuitous due to errors in all of the solutions in terms of accurately calculating the turbulent kinetic energy. 

We then investigated the SWTBLI cases which were the motivation for this work. Experimental data for 
the UFAST test case, with Mach 2.25 flow approaching the oblique shock interaction, indicated that —u'v'/k 
exceeded 0.31 at the beginning of the SWTBLI and further downstream near reattachment. Also interesting 
was that the profile nearest the center of the interaction showed the peak value of —u'v'/k to temporarily 
drop. From the computations made with various values of a\ in the SST limiter, the case with a\ = 0.355 
most closely matched the separation size from the experiments, and this choice for a\ was consistent with 
the turbulent statistics from the experiment. A final case investigating a complex system problem of a direct 
connect scramjet experiment also showed that the SST model with a\ = 0.355 provided the best agreement 
through the key shock-separated flow region. 

As with many settings for key turbulence model coefficients, it is very difficult to find universal settings 
for all of the coefficients. The same is likely the case for the a\ coefficient examined in this study. For the 
incompressible test cases considered, it is difficult to recommend a setting different than a\ = 0.31 that 
is in the original formulation. However, it remains somewhat troubling that for the axisymmetric diffuser 
case, which was one of the showcase problems to recommend usage of the SST model, that the default value 
for a\ is not consistent with experimental observation. For the range of SWTBLI cases investigated in this 
study, utilizing a value for a\ closer to 0.355 yields improved predictions and is consistent with experimental 
observation for mean flow and turbulent statistics. 
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Figure 1. Skin friction for incompressible flat plate case. 
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Figure 2. Velocity profile for incompressible flat plate case at Re x = 4.3e6. 



Figure 3. Eddy viscosity profiles at Re x = 4.2e6. 
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Figure 4. Turbulent shear stress profiles at Re x = 4.2e6. 



Figure 5. Wall pressure comparison for Driver axisymmetric diffuser. 
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Figure 6. Skin friction comparison for Driver axisymmetric diffuser. 


NAS A/TM— 20 13-217851 


13 


y (mm) y (mm) y (mm) y (mm) y (mm) y (mm) 



-300 -200 -100 


0 x (mm) IOC* 200 300 


(a) Experiment 



-300 -200 -100 


°x(mm) 100 200 300 


(b) SST 



-300 -200 -100 


°x(mm) 100 200 300 


(c) SST, a-i = 0.34 



-300 -200 -100 


°x(mm) 100 200 300 


(d) SST, ox = 0.355 



(e) SST, ai = 0.37 



(f) BSL 

Figure 7. Turbulent shear stress contours for Driver axisymmetric diffuser. 
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Figure 8. Behavior of —u'v'/k for Driver axisymmetric diffuser. 
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Figure 9. Turbulent kinetic energy contours for Driver axisymmetric diffuser. 
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Figure 10. Skin friction coefficient for backward facing step. 



(a) SST 



Figure 11. Eddy viscosity contours for backward facing step. 
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Figure 12. Structure parameter for backward facing step. 
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Figure 13. Axial velocity contours for UFAST test case in vicinity of SWTBLI. 
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Figure 14. Axial velocity profiles for UFAST test case in vicinity of SWTBLI. 
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Figure 15. Behavior of —u'v'/k for UFAST case in vicinity of SWTBLI - taken from data of Ref. 29. 
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Figure 16. Skin friction for Mach 5 SWTBLI. 
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(a) Mach number contours along cut plane through fuel injectors 



(b) Pressure comparisons along flowpath centerline 

Figure 17. HIFiRE HDCR scramjet at Mach 5.8 freestream operating point. 
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